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The conversion of neutron star to strange star is argued to be a two step process. The first process involves the 
deconfinement of nuclear to two-flavour quark matter. The GR results shows that the propagating front breaks 
I up into fragments which propagate with different velocities along different directions. The time taken for this 

, conversion to happen is of the order of few ms. This calculation indicates the inadequacy of non-relativistic (NR) 

^ ' or even Special Relativistic (SR) treatments for these cases. 

^ . PACS numbers: 



I. INTRODUCTION 

' ' Strange Quark Matter (SQM), consisting of approximately equal numbers of up (u), down (d) and strange (s) 
2 ' quarks, is conjectured to be the true ground state of strong interaction [U, SQM could naturally occur in the cores 
of compact stars, where central densities of about an order of magnitude higher than the nuclear matter saturation 
I density. The transition from nuclear (hadronic) to quark matter should proceed through a conversion to an initial 
stage of (metastable) two flavour quark matter, which should decay to the stable SQM. Thus, neutron stars with 
sufficiently high central densities ought to get converted to strange, or at least hybrid, stars. These transitions could 
have observable signatures in the form of a jump in the breaking index and gamma ray bursts [1, 
' There are several plausible scenarios where neutron stars could convert to quark stars, through a "seed" of external 
, SQM Q , or triggered by the rise in the central density due to a sudden spin-down in older neutron stars_0. Several 
OO ' authors have studied the conversion of nuclear matter to strange matter under different assumptions 0, S Q • 

In our recent work we have argued that the conversion process is a two step process. The first process involves 
\^ > the deconfinement of nuclear to two-flavour quark matter and the conversion process takes some milliseconds to occur. 
' GR effects in such processes is studied by us [TT|. In this article, we make a detailed study of this process. 
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II. THEORY 



We use the Nonlinear Walecka model for the nuclear matter euqation of state (EOS). We consider the conversion of 
' nuclear matter, consisting of only nucleons {i.e. without hyperons) to a two-flavour quark matter. The final compo- 
sition of the quark matter is determined from the nuclear matter EOS by enforcing the baryon number conservation 
during the conversion process. 

The metric describing the structure of the star, is given by [l3 

ds^ = -e^+^dt^ + e2"(dr2 + r^dO^) + 

e'^'Pr^sin^e{d(t>-ijjdtf (1) 

The four gravitational potentials a, 7, p and u) are functions of 9 and r only. The solution of the star is obtained from 
the 'rns' code (13i] . 

We heuristically assume the existence of a combustive phase transition front of infinitesimal thickness, and study 
the outward propagation of the front. Let us assume that the conversion front is generated at the center of the star. 
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FIG. 1: Variation of velocity of the front along the radial direction of the star for three different cases, namely SR, non-rotating 
GR and rotating GR. 
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FIG. 2: Variation of velocity of the front along the radial direction for different rotational velocity of the star. 



and it propagates outwards through the star with a certain velocity. Employing the conservation conditions [l^ and 
further employing the entropy condition [ISj , we determine the flow velocity of matter in the two phases vi and V2- 

It is possible to classify the various conversion mechanisms by comparing the velocities of the respective phases 
with the corresponding velocities of sound in these phases. For the conversion to be physically possible, velocities 
should satisfy an additional condition, namely, < < 1. If we plot different velocities against baryon number, we 
get curves from which we could calculate the initial velocity of the front at the centre of the star The results 
[lOj also show that the range of values of baryon density, for which the flow velocities are physical, increases with 
temperature. Starting with this initial velocity, we investigate the evolution of the front with time. Treating both 
nuclear and strange matters as ideal fluids, the system is governed, together with the metric and the EOS, by the 
Einstein's equation Ri'^ — ^Si'^R = nTi^ and the equation of motion T^'i-k = 

The above two equations are the starting point for deriving the appropriate continuity and Euler's equations pd| . 
After a bit of algebra, we get a single differential equation for v: 



dv 
Or 



W'^v[K + Kl] 



2[i;2(l + G)2-n(l + t;2G)2 



(2) 



uj = and sin6 = 1 in this equation yield the equation for the static star and if we put all potentials equal to zero, 
we recover our equation for the SR case p^ . 
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FIG. 3: Variation of velocity of the front along the radial direction for different X- 
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FIG. 4: Variation of time of arrival of the front at certain radial distance for different cases. 

III. RESULTS AND DISCUSSION 

Having constructed the density profile of the star for a fixed central density, the respective fiow velocities vi and 
V2 of the nuclear and quark matter in the rest frame of the front , at a radius infinitesimally close to the center of the 
star. This would give us the initial velocity of the front {—vi), at that radius, in the nuclear matter rest frame. With 
this, we integrate eqn. (2) outwards along the radius of the star. The solution gives the variation of the velocity with 
position as a function of the time of arrival of the front. Using this velocity profile, we can calculate the time required 
to convert the whole star using the relation ^ = vG. 

For a rotating star, due to the asymmetry, we introduce a new parameter x = cos6, along the vertical axis of the 
star. We start our calculation by choosing the central density of the star to be 7 times the nuclear matter density, 
for which the Keplerian velocity of the star is 0.67 x 10~'*sec~^ (the rotational velocities given in fig. (2) are all in 
units of 10~^sec~^). For this central density, the initial velocity of the front comes out to be 0.45. In fig. (1) the 
propagation velocity of the front along the radial direction of the star for three cases. The unbroken curve is for the 
SR case, the broken curve for non-rotating GR case and the dotted curve for the rotating GR case with x = 0, i.e at 
the equator. Due to the asymptotic behaviour the velocity shoots up at the centre and saturates at larger radii. It 
can also be clearly seen that the GR effect increases the velocity of the front considerably (max;imum by 30%) and 
the effect is most pronounced for the static case. The rotational effect of the star seems to suppress the GR effect 
and therefore the velocity of the front decreases. The result becomes clearer if we look at fig. (2) where we have 
plotted the front velocity with equatorial radius for different rotational velocities; as the rotational velocity increases, 
the velocity of the front decreases. 

From fig. (3), we find that the front velocity is maximum along the polar direction and minimum along the equator. 
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Therefore, at any particular instant of time, we may have a situation where the polar part of the star has been 
converted while along the equatorial direction, the front is still propagating. 

Prom fig. (4) we can see that the time taken by the conversion front to convert the neutron star to two-flavour 
quark star is of the order of few ms. The static star takes the minimum time (3.3tos) whereas the rotating star takes 
the maximum time (5.1ms) due to the enlarged equatorial radius. The polar part of the star needs much lesser time 
for conversion (3.1ms), even less than static star, as its radius gets compressed. 

To summarize, we have shown in this article that the conversion of nuclear matter to quark matter in compact stars, 
especially rotating stars which are more realistic than static stars, is strongly affected by GR effects. The emergence 
of different conversion fronts, propagating with different velocities along different radial directions was not anticipated 
by Newtonian or SR calculations. It remains to be explored whether the incorporation of dissipative effects materially 
changes the results. Though the calculation is much involved, such an investigation is on our immediate agenda. 
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